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ABSTRACT 


A seismic signal processor has been developed and tested for use vtrith 
the NOAA-OOES satellite data collection system. Performance tests 
on recorded, as well as real time, short period signals indicate that 
the event recognition technique used (formulated by Rex Allen) is 
nearly perfect in its rejection of cultural signals and that data can be 
acquired in many swarm situations with the use of solid state buffer 
memories. Detailed circuit diagrams are provided. The design of a 
complete field data collection platform is discussed and the employ- 
ment of data collection platforms in seismic networks is reviewed. 
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APPLICATIONS OF SATELLITE DATA RELAY 


TO PROBLEMS OF FIELD SEISMOLOGY 


INTRODUCTION 

In 1975, The Geophysics Branch of Goddard Space Flight Center decided that developing a data 
collection platform (DCP) to transmit seianic information by satellite relay would be an excellent 
way of using space techniques to benefit scientific research. For the next several years, a coopera- 
tive program was maintained with Rex Allen at the U.S. Geological Survey’s Branch of Earthquake 
Mechanics and Prediction at Menlo Park, California, with Goddard furnishing partial financial sup- 
port for the development of a seimiic-event detector algorithm. In 1977, a breadboard event 
detector, using a first version of Allen’s algorithm, was designed, built and tested by Rovert Novas 
(1977)* at Goddard. This was preliminary to the present effort. 

The des^ goal was a system with maximum reliability and scientific return at minimum unit cost 
and complexity. Scientific requirements were established by a survey of potential users in univer- 
sities and federal agencies. 

Upon completion of the initial design in mid- 1977, it was decided to construct a breadboard engi- 
neering model to demonstrate the viability of the concept. The purpose of this effort was to show 
that a satellite seismic DCP can be constructed with no technical risk. All of the elements of the 
DCP that might represent a risk were breadboarded and the results were used to refine the final 
design of the DCP. The breadboard was completed in early 1979 and was initially tested using 

*Novu, R. G., 1977, An Application of Microprocessor Technology to Remote Station Analysis of Seismic Signals, 
unpublished Matter of Science Thetis, Lehigh University, Bethlehem, Pennsylvania. 
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magnetic tapes of Alaskan seismic events, ftimished by the University of Alaska. Since then the 
breadboard has been undergoing operational tests using a seismic signal transmitted by telephone 
lines to Goddard from a vertical axis, short-period seismic installation near Baltimore, Maryland. 

The stated goal of proving the feasibility of the concept has been accomplished. Programmatic 
considerations have precluded further efforts to use the existing unit or develop a field-hardened 
unit. The purpose of this report is to describe and evaluate the breadboard design and operational 
characteristics. Additional information may be furnished to any group desirous of continuing this 
development. 

HISTORY AND JUSTIFICATION 

The collection of data by satellite is a relatively new technique Hrst demonstrated in 1967 using 
NASA's ATS-1 (Applications Technology Satellite). The first demonstration was the NASA Omega 
Position Location Equipment System which proved that accurate positions could be obtained from 
platforms in remote locations and that satellite relay did not degrade the data. This experiment was 
followed in 1969 by the Interrogation, Recording and Location System flown on Nimbus-3 and 
Nimbus4. This was the first global satellite system to demonstrate the worldwide capabilities of 
data collection by satellite. 

Because they were designed to respond to interrogations from the satellites, these ground systems 
were relatively large and expensive and required considerable power. This was overcome in the 
Landsat series of satellites, initiated in 1972, by designing the ground platforms to transmit at ran- 
dom tunes, thus eliminating the requirement for having a receiving system in the DCP. In 1974, 
NOAA introduced the (30ES (Geostationary Operational Environmental Satellite) system that 
employs either a scheduled or satellite interrogated transmission system. 
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Figure 1 it a block diagram of typical latellite relaying systems currently being used to return infor* 
mation from low<data*rate geophysical instruments such as tide gauges, strain meters and tiitmeters. 
However, because of the high*data-rate requirements, no practical cost-effective system presently 
exists for returning high-data-rate seismic information. For example, continuously recording seis- 
mic data, using a 12-bit word for signal resolution and sampling at SO hertz (Hz), requires almost 52 
megabits per day per component. 

Off-the-shelf availability of such a field system would have many advantages. Currently, most un- 
manned seismic field stations either have to be visited every day or two, to replace recording paper, 
or the information has to be transmitted to a central location by expensive and sometimes noisy 
phone lines and/or radio relays. Phone lines, almost nonexistent in remote, inhospitable or under- 
developed areas mch as Alaska, are often unreliable, even in populated areas. Furthermore, ground 
communcations often become inoperative before, during, and after a major earthquake. When geo- 
physical systems are operated in extremely inaccessible regions, data are usually preserved on low- 
powered, slow-speed recording systems which may run unattended for months; the data are then 
collected several times a year. Such systems require sacrifices in timing accuracy and information 
content, and since data analysis is often delayed for months after the events, earthquake prediction 
capability is lost. Also, there can be no assurance that the instrument is performing as planned. In 
addition, it is often desirable to rapidly augment a seismic network to collect earthquake precursor 
signals or monitor after^ocks, and the dependence upon phone lines or radio relays might impede 
the mobility of instrument siting and increase installation time. Finally, for earthquake disaster 
relief, it would be of inestimable value to have available worldwide seismic data in real time. The 
primary disadvantages are the increased cost and complexity of the collection system and. depending 
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Figure 1 . Block diagram of a typical satellite data relay system 



on the requirementi of the inveitigator, the possible necessity of working with simplified or de* 
graded data. 

DESIGN PHILOSOPHY AND USER REQUIREMENTS 

To be most widely applicable, a seismic DG^ should possess the following characteristics: 

1 . Provide, in near real time, significant sc’entific data for a broad spectrum of investigators 

2. Have a reasonable price; i.e., be affordable by most investigators 

3. Operate with existing satellites 

4. Operate on a one-to^ne basis with a single seismic system ; i.e., not be dependent on cross* 
correlation schemes between multiple systems 

5. Be battery operated with at least a six*months life between battery changes 

6. Be field hardened; i.e., reliable, capable of unattended operation, environmentally sealed, 
wide thermal operating range, minimal moving parts such as tape recorders, etc. 

7. Be relatively small and li^tweight. 

The obvious key to a practical field system is an “event detector” device that reliably differentiates 
signals from background noise. Once this is done, the noise periods can be discarded and the event 
signals can be operated on by the system. If desirable, further data compression can be performed 
on the stored events before transmission, igure 2 schematically illustrates such a system. 

The majority of event detecting devices have generally depended on a manually-set threshold for 
comparing short-term energy (signal) with long-term energy (noise). The reliability of such a device 
is considerably increased when cron-correlation between multiple seismic stations is possible, but 
such correlation is obviously not feasible when a single seismometer/DCP system is under 
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considention. MIenby et al. (1977), detailed the development of seismic event detectors. The 
algorithm used for Goddard's breadboard was developed by Rex Allen (1978) of the U.S. Geologi- 
cal Survey in Menlo Park, California, and was based on an earlier design by Stewart (1977). 

In considering the scientific data provided by the system, it was decided that the DCP should be 
applicable mainly to research presently associated with remote, untended, short-period seismic 
installations. This would generally restrict the use of the system to local and regional data studies 
and would avoid data requirements associated with relatively complex analyses of very distant 
events or surface waves. The DCP should then be useful for studies related to: 

1. Crustal Hazard Reductions 

a. Earthquake mechanisms 

b. Earthquake prediction 

c. Interplate and intraplate stress and tectonism 

d. Volcanic eruption prediction 

e. Seismicity of reservoir filling. 

2. Crustal and Mantle Composition and Structure 

3. Mine and Quarry Blast Monitoring 

4. Tsunami Prediction. 

The next scientific design consideration was what components of the individual seismic signal are 
needed for ihe various studies. In order of increasing data complexity these arc; 

1 . Number of events per day Volcano monitoring 

Earthquake swarm studies 
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2. “P” (compressional wave) arrival 

time 

Direction first motion 
Duration and/or maximum amplitude 
and frequency 

3. AU of “2” plus “S” (shear wave) 

arrival time 


Location and magnitude of earthquakes 
(tectonic and volcanic) 

Tsunami prediction 
Blast monitoring 
Fault plane solutions 
Earthquake prediction 
(Vp/V^ anomalies) 

Regional seismicity 


The challenge, then, was to design a practical field DCP system that would provide as much as 
possible of the above information. To help us in this, university and government scientists were 
consulted regarding their data preferences. Initially, some consideration was given to processing 
the data in the field and relaying back only numbers represent'i.iive of the values of times of the 
desired features. I^owever, developing an algorithm to identify the “S” phase would be a ve»y 
difficult, if not impossible, task. In addition, we found almost no application in which the users 
were willing to accept the loss of the actual trace, primarily because of a natural unwillingness to 
depend on a field computer to analyze the signal. For these reasons it was decided to reconstruct 
the returned signal into analog form. The general requirements for such a signal were a bandwidth 
of O.S to 25 Hz, a maximum-event length of around 180 seconds, and a digital resolution of 12 bits 
(72-dB signal-to-noke ratio). Considerable interest was also expressed in using 16-bit word lengths 
for signal level (96-dB signal-to-noise ratio), but, at that time 16-bit analog-to-digital converters 
lacked stability and reliability. For these reasons the final system was designed for a 12-bit word at 
a sampling rate of SO times per second. Thus, a 180-second recording involves a total of 108 kilobit 
(kb) (not including any overhead due to housekeeping, timing, magnitude and quality data). 
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In addition to th^se primaiy requirements, other factors arose. First, because the S arrival from 
nearby events is often stronger than the P on short-period vertical sensors, it is desirable to return 
a portion of the trace preceding the selected event to verify that the event was picked on the 
P and not the S phase. For regional earthquakes, the headwave P^ is weak in comparison with P 
and P*. It is therefore desirable to have, perhaps, 5 seconds of pre-event detect signal so that those 
phases can be properly identified. It would then be possible to put tight constraints on depth and 
distance and provide detailed information on regional structure. 

A short pre-event strip is also useful for indicating the background noise level and hence the opera- 
ting reliability of the event detector. For these reasons, a 10 second pre-event strip precedes the 
recording of the actual event. As mentioned previously, the maximum total recording time per 
event is 180 seconds (including the 10 seconds pre-event time). However, this total time is adjus- 
table because, for many applications, an event time of 90 seconds is sufficient. 

Several users expressed concern about the possible saturation of the system in the event of swarms. 
A number of schemes were considered. A procedure of buffer swapping to be described below was 
adopted. Three such buffers or memories would allow an efficiency of 43 percent in the event of 
swarms. The system would saturate the available output data stream but would be able to record 
43 percent of the time for transmission. 

Magnetic tape was eliminated for event storage because of mechanical complexities. Reliable 
bubble memories are not yet available, and power requirements are high. Solid-state memories 
proved to be quite satisfactory. The breadboard contains two memories. When an event is identi- 
fied, number one memory records 1 0 seconds of pre-event noise and, depending on the setting. 
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80 seconds of the event at a high-data rate. The stored event is then “dumped” at a lower-data rate 
through the satellite. It requires about 9 minutes for a ll^minute event to be transmitted to the 
satellite. 

The design was dependent on the choice of satellites. A dedicated channel on a synchronous satel- 
lite would permit continuous transmission (depending on the power budget of the DCP). In con- 
trast, nonsynchronous satellites require satellite callup, random or timed data dumps. In both cases 
maximum data rates vary depending on antenna sizes, power, etc. While there are numerous com- 
munication satellites that are technically suitable, ground unit costs are related to the operating 
frequencies of the satellite. Thus, our requirement for low DCP unit costs eliminated many satel- 
lites from contention at the present. Most of the high volume satellites operate in the 1-GHz (gig- 
ahertz) or S-GHz satellite allocations. Technology is not yet up to producing inexpensive and 
efficient transmitters at these frequencies. A typical 20-watt transmitter at 2 GHz is about 10 
times as expensive and half as efficient as its 400-MHz (megahertz) counterpart. In addition, 
because frequency slots are assigned within satellite transponders to a high percentage accuracy, the 
frequency control is much more expensive at the microwave frequency than at uhf (ultra-hi^ fre- 
quency). 

Accordingly, satellites with uplink frequencies in the uhf range are preferred. As an example of the 
maturity of the uhf technology, a single-module power amplifier capable of generating a IS-watt 
output signal from 1 SO mW (milliwatt) of drive at 400 MHz costs about S80. A similar microwave 
power amplifier costs $10(X) and is half as efficient. 

The most extensive network of satellites using a uhf data collection system (DCS) is the GOES 
system. In addition to its prime function as an imaging meteorological satellite, GOES has a 
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400-MHz uplink DCS. A simplifled block diagram of the DCS system is given in Figure 3. Note 
that there are 200 DCP uplink channels between 401 .2 and 401 .7 MHz. Each of these channels is 
1 S-kHz (kilohertz) wide and is intended to accommodate ASCII code at 100 bps (bits per second). 
The satellite, being synchronous, allows random dumping by the DCP whenever an event is identi- 
fled and stored. Also, since the United Nations* World Meteorological Organization (WMO) protocol 
provides for a worldwide GOES system, it seems likely that a GOES type DCS will be available for 
at least the near future. 

Therefore, the design and demonstration work was conducted, assuming the GOES DCS character- 
istics as the design driver. However, because the 100 bps restriction is relatively severe (and, in fact, 
represents a “worst-case" situation for all practical purposes), and since microwave transmitter 
technology is fast becoming mature, a modular approach was adopted which would allow an easy 
change of output data rate and transmitter frequency. 

BASIC OPERATING PRINCIPLES OF SYSTEM 

The output of a single-axis, high frequency (1 to 2 Hz) seismometer is continuously monitored by 
an automatic event detector (“P” picker). When an event is identified, up to 180 seconds of signal 
(10 seconds pre-event noise and 1 70 seconds of event) is recorded and stored in a solid-state 
memory at a sampling rate of 50 times per second and a 1 2-bit word for signal resolution. A delay 
line allows the system to recover 10 seconds of pre-event signal after the event picker decides it has 
an event. 

Upon completion of recording, the first memory system goes off the line and begins transmitting to 
the GOES satellite at 100 bps. During the 1 8 minutes required for the first memory to dump a 
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Figure 3. GOES data collection system data flow. 



3-minute signal to the satellite, a second memory is on the line to record the next detected event. 
For fleld use, particularly if swarms are expected, at least three memories would be required. 

DEMONSTRATION HARDWARE 

The demonstration breadboard was designed and constructed according to the following criteria: 


Input signal 

Analog 

Bandwidth 

25 Hz 

Event length 

180 seconds (maximum) 

Resolution 

12 bits 

Output signal 

Compatible with GOES (100 bps bi-phase) 

Operating mode 

GOES emergency event triggered 

Power 

Battery pack 

Battery life 

Six months (average 1 2 events per day) 

Cost per fleld unit 

Less than $10,000, including DCP and radio set 


A block diagram is given in Figure 4 and contains all the subunits of the breadboard. In what 
follows, the design of the demonstration unit on a subunit basis will be discussed. 

The breadboard receives signals from an event simulator, a tape recorder, or a conventional discrim- 
inator. The prerecorded analog tapes were provided by the University of Alaska. The event simula- 
tor generates a damped harmonic signal electronically. The breadboard output is a serial-digital 
signal at the GOES rate of 100 bps. This signal is passed to a digital-to-analog converter for com- 
parison with the input signal. 
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Figure 4 . Block diagram of seismic DCP showing portion cons 





ANALOG CIRCUITS 


Figure 5 shows the analog circuit block diagram, a detailed circuit diagram is shown in Sheet 2 *. 
The instrumentation preamplifier has a bandwidth of SO Hz and a gain of two. This amplifler 
gain can be increased to 2000 by a component change. The high gain was not required for the 
breadboard because the tape recorded signal was already preamplifled. The wide bandwidth of 
SO Hz enables the event detector to determine event-occurrence time to within 10 milliseconds. 
With a SO-Hz infonnation bandwidth, the Nyquist sampling theorem dictates at least a 100 sample- 
per-second rate. The 1 2-bit analog/digital (A/D) converter has additional filtering (2S-Hz low pass) 
to minimize signal aliasing. 

The demonstration breadboard could have been designed with one 12-bit A/D converter followed 
by a digital filter and a divide-by-16 circuit. This approach was not used oecause: (1) using the 
digital divider and filter would have required more modules, and (2) this approach also allows an 
easy change of microprocessors since the entire event detector is isolated from the main data 
stream. 

Digital Delay 

A delay is required before buffer storage to: 

1. Provide the experimenter with some pre-event noise for signal analysis. 

2. Provide time for the microprocessor to calculate whether an event has occurred. 

3. Provide pre-event time for the base station receiving system to obtain synchronization. 

4. Provide time for the DCP transmitter to stabilize prior to sending an event signal. 


*There ue 1 1 engineeting tduepiinti lefened to u theett. (See beck cover.) 
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A delay of 10 seconds appeals to be adequate to perfonn these functions. A delay of 10 seconds is 
obtained by a 1024 word, 1 2-bit*per-word CMOS RAM used as a first in*first out memory. If more 
delay is necessary, the delay time can be dianged to a maximum 8 1 .92 seconds by a wiring modifi- 
cation that changes the decoder input signals (Sheet 3, module £3). The decoded output resets the 
12-stage ripple counter module El at the required time. The output of the ripple counter also 
provides the addressing signals to RAM devices with address-state changes every 20 milliseconds. 

Data delayed by the delay time is always being sent to the buffer memory module for possible 
storage. Recording in the buffer storage depends on the buffer storage control and gating signals 
that are under control of the microcomputer. 

100 kb-Buffer Mwnory 

There are several solid state technologies that can be used for a memory size of 100 kb; i.e., core, 
plated wire, CMOS, NMOS, MNOS and magnetic bubble. Because of reliability considerations, 
mechanical devices (i.e., tape recorders) were not considered. Charge coupled devices (CCD) were 
not a candidate because of our need for a low operating rate of 100 bps. The operating rate has to 
be greater than 50 kbps for most (XD chips containing 4096 bits or greater because of “dark cur- 
rent" limitations. Plated wire is too expensive; core dissipates too much power compared to the 
other solid-state devices. Metal-nitrate-oxide semiconductor (MNOS) devices are too expensive, 
ease of manufacture is poor, and the availability of second sources is also poor. A list of the candi- 
date components and their characteristics for a 100-kb buffer memory is given in Table 1 . 

The static CMOS RAM memory device was selected over the other candidates primarily because of 
the very low average power dissipation. Rejection of the magnetic bubble device was not primarily 


Table 1 

Memoty G>mponentt Characteristics for Mass 
Memory Application 



CMOS 

DYNAMIC 

NMOS 

MAGNETIC 

BUBBLE 

STATIC 

NMOS 

Manufacturer’s Number 

1M6S03 

Intel 21 16 

TBMOlOO 

EMM 4044 

Chip Density (bits) 

1024 

16384 

92000 

4096 

Chip Organization 

RAM 

RAM 

FIFO 

RAM 

Number of Different Voltages 
Required 

1 

3 

4 

1 

Operating Temperature (*C) 

*55+125 

*55+125 

+15+35 

*55+125 

Module Cap., kilobits 

110.6 

114.7 

92 

110.6 

Average Power, watts 

0.054 

1.7 

1.1 

6.4 

Module Component Cost (S) 

540 

390 

470 

416 


due to power requirements. The magnetic bubble device has the advantage that storage is nonvola* 
tile and the device can be power switched. We did not use the device because its availability is poor 
at present and the ease of use in the design is difficult. A re-evaluation of the magnetic bubble 
should be conducted in a few years, when the device’s performance, availability, and adaptability 
are improved. 

A detailed schematic for each of the two 1 10*kb buffers is given in Sheet 7. Note, that the buffer 
is parallel organized (12 bits X 9216 words). The storage elements are CMOS RAMS (1024 X 1). 
Operating length of each buffer can be changed by eiglit switches from 2048 to 9216 words in steps 
of 1024 words. This corresponds to event lengths of 41 to 184 seconds. (S-P range circles of about 
250 kilometers (km) to 1 700 km.) 


18 



E^ftt Dttictor 


The event detection function it peifonned by a microprocenor that it programmed to proce^ 
digitized teiimic tignali in real time. Interface to the rest of the system is particularly simple. 
Figure 6 shows the event detector interface signals. The only output interface signal used by the 
remaining modules it the event-detect signal; the event*status data, which would be used in a field 
unit, it not used here; this signal can be obtained via the data terminal. The control signals and 
program constants are on front*pa>’el switches and are read only during program initialization. 

There are several microprocessor systems that could have been used. A list of candidate micro- 
processors is given in Table 2. The CMOS CDP*. 802 was selected based on the following criteria; 

1 . I-">w power 

2. Add time 

3. Support chips available 

4. Reliability (only microprocessor on GSFC preferred parts list). 

The CDP1802 microprocessor is a single-chip, 8-bit, static microprocessor fabricated in CMOS 
technology. The CDPl 802 thus hf i all the advantages of CMOS technology ; i.e., low power dis- 
sip?.<-=>n, sin<;ii.' wide-range power supply, full operating temperature range and a single-phase clock. 
Our system uses a i-volt powrr supply and a 2-MHz clock. (With a 2-MHz clock, the machine-cycle 
time is 4 mi.roscconds and the instruction cycle time is 8 microseconds.) 

Refer to Sheets 4, S, and 6 for the event detector module details. Note that Sheet 4 shows the 
interface and control circuits and the connections to the CDPl 802 microprocessor. Sheet 5 dia- 
grams the 4K words of RAM (used as working storage) (32, CDPl 822’s). Sheet 6 shows the 4K 
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PORTABLE 
DATA TERMINAL 




Table 2 

Candidate Microprocetion Characterization 
With a 2*MHz aock 



Access 

Power 

(Milliwatts) 

Data Bus Width 

Add Time 
(Microseconds) 

SBP 9900 

1*L 

500 

16 

9 

TMS 9900 

NMOS 

1000 

16 

IS 

CDP1802 

CMOS 

4 

8 

8 

1 M 6100 

CMOS 

5 

12 

10 


words of ROM (4, M 2700) and 32 words of RAM used for program and register storage). Only 2K 
words of ROM are actually needed ; l K words for a standard utility program and 1 K words for the 
seismic processor program. The utility program is used to communicate with the data terminal. 

Control Circuitry 

The control circuitry is used to generate; 

1 . Analog-to-digital conversion sample pulses 

2. Buffer<ontrol signals 

3. Multiplexer-control signals. 

All of these signals are derived from combinations of the 2-MHz clock; the event-detect signal and 
the buffer full/empty signal. 

The 2-MHz clock is divided down to generate all the sampling and timing pulses (modules Cl . CIO. 
Sheet 8). Also derived are the analog-to-digital converter (ADC) sampling pulses which are contin- 
uous at either SO cps (1 2-bit ADC) or 100 cps (8-bit ADC). 
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Figure 7 shows the buffer timing sequences that are generated on the control circuit board. Until 
an event-detect signal occurs, both buffers are in standby. When the event detector declares a valid 
event, the event-detect signal triggers buffer 1 into operation. Buffer 2 remains in standby. Soon 
after the event-detect signal occurs, the buffer- 1 clock starts operating at its high rate (SO cps) and 
the MWR-1 signal enables a write operation. Buffer-1 initialization occurs when the first clock pulse 
causes memory location one to be written into. After buffer 1 is full, the “1 full” signal is gener- 
ated. The buffer-full signal starts a read operation clocked at the GOES rate (100 bps). This is 
done through the MWR signal that places the memory into a read state. The memory addressing 
is organized so that the first-clock pulse after full signal enables reading from memory location one 
(the memory is a fiist-in/first-out type). After all memory cells are read, a (1 empty) signal is 
generated which places buffer 1 into standby. However, if another event occurs between the buffer- 
full and the buffer-empty signal, buffer 2 begins a write operation. Buffer 2 will not perform a read 
operation until buffer 1 has received an empty signal, and a buffer-2 full signal is generated. 

The data-ready signals, diagramed in Figure 7, control the multiplexing of the two buffer output 
signals into one signal during buffer read times (Sheet 9). In a field DCP, this signal would be bi- 
phase modulated and then sent to the transmitter. In our unit, this signal is sent to a digital-to- 
analog converter, then to a visual recorder. The recorder used is a standard Sprengnether three 
channel drum recorder. The three signals recorded during unit-performance testing are: (1) the 
analog signal after the preamplifer, (2) the event-detect signal, and (3) the delayed processed data 
from the buffers read out at the equivalent of 100 bps. 

The unit was constructed of CMOS DIP integrated circuits mounted on Augat circuit boards with 
connections by wirewrap. These boards are mounted into a standard 48-cm (19-inch) rack chassis. 


22 



23 


Figure?. Buffer timing diagram. 


8.S cm (3Vi inches) high. Power supplies are mounted separately. The unit was partitioned into 
boards as follows: 

Board Function 

1 Analog circuitry, control logic and multiplexing 

2 Microprocessor, associated random access memory and 

read only memory 

3 Buffer memory # 1 

4 Buffer memory #2 

The system was partitioned so that: (1 ) a different type microprocessor module could easily be 
added for additional evaluation, and (2) buffer memory could easily be expanded if necessary. All 
controls are front-panel mounted with exception of the buffer-length switch. 

EVENT-DETECTOR ALGORITHM AND ITS IMPLEMENTION 

The event-detection program, used on the 1802, is based on an algorithm developed by Rex Allen 
(1978) for the automatic detection and timing of seismic events from a single seismometer; how- 
ever, modifications were necessary to run the program on the 8-bit RCA 1 802 microprocessor. The 
program is an interrupt-driven (real-time) task that identifies events to within 10 milliseconds. The 
program also evaluates the accuracy of its picks, thus eliminating the recording of events generated 
from noise sources such as vehicle traffic. 

Appendix 3 contains the 1802 assembly code, a memory map and the tables for conversion of con- 
trol constant values to switch settings. Appendix 2 contains a running description of the Allen 
(1978) algorithm as implemented for the 1802. 
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Data from the 8-bit analog/digital converter is searched for the possibility of an event according 
to Allen’s criteria. The characteristic-function calculation is the primary time consumer of this 
event-search mode operation. Ideally, the whole event-search process for one sample should not 
take more than 10 milliseconds. In practice, the average time was calculated as 9.64 milliseconds, 
and in the worst case, 16.4 milliseconds. In actual use, we found that most samples (85 percent) 
could be handled in 10 milliseconds. When the program requires more than 10 milliseconds to 
process the sample, the next sample is ignored. The consequences of this time constraint are dis- 
cussed in the en^eering tests section. 

Once a potential event is registered, the program enters the event-validation mode to test whether 
the suspected event passes duration, frequency, and amplitude criteria. On the average, this process 
should take 6.52 minutes, and in the worst case 13.28 milliseconds. In practice, we have not 
observed undersampling during the event-validation mode. 

The current formulation of the algorithm will store up to 256 event initiation times in the form of 
clock cycles since initialization. Interpolation to a fractional clock cycle is not done. In addition 
to the event times, the zero crossings and peak amplitudes used in the analysis (up to 1 28) are also 
available. The memory lap in Appendix 3 shows that of the available 4K words of RAM, only 
about 1 .5K words are used. Of the 4K words of PROM, IK words are required for the event- 
detection program, while IK words are used for the utility. This small RAM/ROM requirement 
indicates a possibility of sharing the memory resources between two processors to decrease the 
apparent cycle time. Although this alternative might permit faster processing of individual samples, 


multiprocessors have not been explored here. 


The expected time requirements during each of the processor modes could be improved by using 
high-speed multiply /divide chips as peripheral devices to the 1802 processor. Althougli 1802 
processor compatible forms of these chips are not yet available, it seems likely that such devices 
will appear in the immediate future. 

Engineering Evaluation 

The completed unit was subjected to several tests designed to evaluate the system’s ability to 
detect events over a wide range of input-signal parameters. These measurements were then used to 
calculate: (1 ) the probability of false detection on broadband noise, (2) the probability of under 
sampling and, (3) the time to detect and verify an event. 

The sensitivity to noise was measured with a “white” gaussian noise signal (i.e., no impulse noise). 
The probability of a false detection and the ability to complete the search and validate tasks in the 
prescribed time were measured. The statistical behavior of this type of noise exceeds, within the 
design bandwidth, the current implementation of the program’s validation test (20 zero crossings 
in 2 seconds). Theoretically, additive gaussian noise in a 50-Hz bandwidth should have 38 maxima 
per second and 50 zero crossings. 

A graph of the probability of false detection as a function of noise and gain is given as Figure 8. 
The noise levels are in millivolts per square root hertz measured after the preamplifier (b»:iuwidth 
SO Hz). Also plotted along the abscissa are the aralog/digital converter quantization levels. With 
the program gain set to maximum, the figure shows that the probability of a false trigger increases 
significantly as the noise level rises above one quantization interval. For this reason, the noise level 
should be adjusted to less than S millivolts. 
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PROBABILITY OF FALSE EVENT DETECT 



Probability of false detection as a function of noise arKi CPU gain. 










The predicted number of instructions for the program to process one sample in the search and 
validate modes, with a full scale signal of 1 27, is as follows: 


Mode 

Worst Case 

Average 

At Best 

Search 

2050 

1205 

795 

Validate 

1650 

815 

705 


For a 2>MHz clocking rate (which yields 8 microseconds per instruction time) there will be 1250 
instructions within the 10 millisecond data sampling interval. It appears that under worse case 
conditions the central processing unit (CPU) data input could be under sampled, since the max- 
imum number of instructions per sampling interval is 1250. Measurements were made to check the 
efficiency on real-time data. The algorithm efficiency was measured by counting the CPU external 
flags (EF) that request transfer of data into the CPU. EF pulses were counted over a 10-second 
period for various input noise and signal levels. Over a 10-second period there should be 1000 
transfers. 

The following results were obtaine J: 

1. No under sampling was detected with noise levels of 5, 10, 20, and 30 mUlivolts. With the 
CPU programmed for maximum gain, the number of instructions to process these noise 
levels varied from 2 to 75 per EF sample time. The number of instructions was directly 
related to the noise level. As expected, the maximum gain setting yielded the largest run- 
ning time. 

2. No under sampling was detected when a simulated event was used as a signal. The simula- 
ted event is an electronically generated damped harmonic sinewave with a natural frequency 
of 10 Hz and a decay time of about 3.0 seconds (Figure 9). Under sampling was checked 
with various peak signal levels (4, 2 and 1 volt) and the CPU gain. 
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(a) SINGLE EVENT 








(b) COMPOUND EVENT 


Figure 9. Electronically generated event signals, vertical axis 2.0 v/cm, 
horizontal axis 0.5 sec/cm; a. single arrival; b. double arrival (P,S). 






3. Under sampling wu detected when a compound lignal was used. A compound signal is 
two simulated events that occur within one second of each other ;i.e., double arrivals or P 
and S phases (see Figure 9b). Under sampling was detected when the CPU was operating 
at maximum gain and the peak>input signal was 4 volts. Under these conditions the prob- 
ability of missing a data transfer into the CPU was calculated to be 0.01 1 . 

4. Figure 10 shows the relationship between the start of an event and the event-detect signal. 
Measurements obtained from over 200 trials, with both simulated and compound signals, 
show a variation from 0.1 to 3.0 seconds from the start of an event to the leading edge of 
the event-detect signal, t^ . Figure 1 1 is a plot of the mean time to detect and verify a 
simulated signal t^ , as a function of peak-signal level and the CPU gain. These tests showed 
that a 10 second delay is sufficient to provide the computational time needed to determine 
that a noncultural event has occurred. 

PERFORMANCE TESTS 

As part of the evaluation, the processor unit was connected to a realtime analog signal from an 
auxilliary short-pt. .nod vertical seismometer (SPZ) at the Geophyacs Branch, Ellicott City, Maryland, 
seismic station (identification code ECM). This signal is unflltered and has a bandpass appropriate 
to the seismometer and voltage-controlled oscillator responses (O.I to 30 Hz). The ECM station is 
located just east of the intersection of 1-40 and 1-70 west of Baltimore, Maryland. The instrument 
package is on a poured concrete slab in contact with the banded member of the Baltimore Gneiss. 
Because of its location, unflltered signals from ECM contain large numbers of heavy and light 
vehicle signatures. 1-70 is, in fact, one of the main truck routes into Baltimore. Local earthquakes 
and close mine blasts are not frequent enough to provide the kind of detailed evaluation of the event 
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Figure 10. Timing relationship between the start of a simulated 
event (a) and the event-detect signal (b). Vertical 
scale 5 v/cm, horizontal scale 0.5 s/cm. 
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Figurt 11. Mean time to detect an event at a function of signal level 
and CPU gain. 







reception aliorithm on wisnic tignalt as reported by Alien (1978). The relatively low frequency 
of luch events only allowt us to show that the Allen algorithm performs on such blasts and earth- 
quakes in a manner consistent with Allen’s observations in a much more active seismic environ- 
ment. An example of the performance on noncultural signals is given in Fi^re 1 2. 

Within a typied seven day period, four mine blasts closer than 200 km are detected at ECM. Also, 
there are two known areas of low-level seismicity within 600 km of ECM (see Bollinger, 1973 and 
Sbar uid Sykes, 1973, for example). The Lancaster, Pennsylvania, area has been responsible for 
about three events per year and the magnitudes are typically well under mb * 3.0. The second 
zone is in west central Vir^ia and has averaged two events per year, also with magnitudes well 
under 3.0. 

Although the low number of close noncultural signals does not permit the quick acquitition of 
meaningful performance statistics on such sisals, the extensive vehicular traffic near ECM allows 
a definitive analysts of the Allen algorithm’s ability to discriminate between cultural signals and 
real seismic events. To illustrate: the unfiltered SPZ output contains about 40 truck and 20 
automobile signatures in a typical 4-hour period centered around one of the “rudi hours." There 
is also a clear diurnal cycle in the high-frequency background noise with peaks during the rush 
hours. As would be expected, there is a strong decrease in the high-frequency background and in 
the number of high-frequency signatures on weekends and holidays. 

Following is a summary of the performance of the system on noncultural signals. The results 
are presented as the probability that the algorithm will fail to detect an event in a particular 


EXAMPLES OF SIGNAL PROCESSING 
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FiMir*12. ExMiiptw of signal proconnsfavtlMbrwdbiMnl unit. Tha 
corr«N>onding to tha avants di^layad art marliad wHh trii 


ditUAce nnfe are baaed on a total of 80 eventi. 


• False hita/poasible false hits 1 0 to 25 percent 

• Activity level of 100 bps channd 20 to 40 percent 

• Event miss probability D>600km 0 percent 

1000 km >D> 600 km 25 percent 
D> 1000 km 100 percent 

In our formulations, Allen’s (1978) algorithm has a soft cutoff at around 1000 km due to our 
choice of averaging technique and the accentuation of the high frequency sensitivity by the form 
of Allen’s characteristic function. Since the characteristic frequency of seismic signals decreases 
with increasing distance, the algorithm will not respond to distant signals unless they are stronger 
than normal. 

The best indication of the performance of the system on cultural signals is to measure the percent- 
age of time that the 100 bps output is active. In general, the activity level is a measure of the rate 
of false triggering. Also, the truck signature is of characteristic frequency greater than 10 Hz and 
shows two sfiort duration q>ikes of about twice the amplitude of the rest of the **<mature. T ’ ^ 
signature provides all the necessary elements for a severe test of the rejection of cultural s*; nals and 
occurs often enough to yield good statistics in a reasonable time. 

Activity on the 100 bps channel is, of course, heavily dependent on the choice of opera; >iig con- 
stants input to the microcomputer. With an optimum set of con^^ants, the activity le* . i 's between 
10 and 20 percent. The same set of constants was used to generate the performance statistics on 
noncultural signals reported above. Since the noise environment was intentionally rrtjde more 
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severe (no filtering) than would normaUy be the case, the rejection of cultural signals should be 
nearly perfect in most applications. 

Since the number of noncultural signals is so small, we tested the buffer swapping procedure with 
cultural signals. A set of constants was selected to give a pick on each vehicle signature with a 
signal-to^ioise ratio greater than a >. The buffer lengths were set for 90 seconds and the system 
was operated through a complete lush-iiour peak (ISO minutes). The number of vehicle signatures 
with suitable signal-to-noise ratio was compared to the number of event-detect signals and the 
number of vehicle signatures in the 100 bps output. Since there are about five signatures per 
minute at the peak of the rush hour, the system was operating at the saturation level for two buf- 
fers. The predicted efficiency is 28 percent and the observed efficiency was 29 percent. The differ- 
ence is probably due to the nonrandomness of the time of occurrence of the signatures. 

Desiffli of a Complete Seismic DCP 

Figure 13 shows five remote data collection platforms and a central data collection station. This is 
a basic form of a GOES-based seismic data collection system. Each DCP is event triggered and uses 
a single DCS channel. The DCP radio sets are small 402-MHz transmitters which have a signal 
bandwidth corresponding to 1 00 bits per second. 

The DCP required EIRP is 48 dBm to communicate at a bit-error rate of 10 ® . A 10-watt trans- 
mitter with an antenna gain of 8 dB is adequate. Figure 14 shows the central station’s received 
signal processing line. This is a low technical risk area since there is nothing unique at the receiver 
as all components have been proven under operating conditions. 


36 



Figure 13. System diagram showing a five OCP network with a central station. 




Received signal processing ai 







Figure 4 is a block diagram of the full DCP. The event-detector output signals to be telemetered 
along with event waveform are the event-detect signal, direction of first motion, event confldence 
measure, and number of timing pulses from event first zero crossing until the event-detect signal 
acknowledges an event has occurred. The DCP transmitter will be power switched under control 
of the microprocessor to increase the battery life time. 

A complete data collection platform unit would consist of: 

1. Seismometer 

2. Event detector with dual buffering 

3. Transmitter (GOES compatible, 10 watts) 

4. Cross-yagi antenna 

5. Battery power pack. 

A power profile was calculated for six months operation. This profile assumes an average of 1 2 
events per day with the DCP/transmitter operating for 1 260 seconds per event ( 1 80 seconds 
record and 1080 seconds playback). Also included in the power calculation are the following 
voltage and current requirements: 

Transmitter (10 watts) 1 2.5 volts at 5 mA (milliamp“-re) idle, 2.5 A trans. 

Event detector 1 2.0 volts at 5 mA 

5.0 volts at 12 mA 
-5.0 volts at 2 mA. 

Power profile dictates that the battery power pack should have the following capacity : 

1 2 volts at 1 900 ampere-hours 


5 volts at 430 ampere-hours. 


It would be posable to derive the $ volts from a 12-volt pack. The estimated cost for a DCP, ex- 
cluding the seismometer, is: 

Event detector and dual buffers $ 1 3,500 

Transmitter and antenna 3,300 

Battery pack 500 

Total $17,300. 

The cost is based on using: 

! . Ceramic CMOS integrated circuit packages 

2. Dual 108-kilobit buffers 

3. Single unit cost (i.e., no quantity discount) 

4. Wirewrap construction. 

The cost of the processor could decrease by as much as $6,000 in large quantities. This cost savings 
would appear as lower costs for the CMOS parts, testing, and packaging. Also, printed circuit pack- 
aging techniques could be used instead of the more costly wirewrap boards. 

The cost of the transmitter and antenna assumes using a HANDAR 524A SMS/GOES data collec- 
tion platform and a high gain crossed yagi. The HANDAR unit contains a GOES compatible for- 
matter, 10-watt transmitter, and power conditioning circuitry. The cost shown above is for a single 
unit procurement. For a large procurement (greater than 10 units) the total cost of a DCP should 
drop to $11,000. 
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OPERATION OF AN EVENT LOCATION NETWORK 
USING THE SEISMIC DCPs 


To illustrate the use of the seismic DCP, we examine the implementation of a location network 
for large events in South America. The purpose of this section is to show how a small network of 
DCPs will allow the location of potentially damaging earthquakes with sufficient dispatch to allow 
mobilization of civil disaster forces should the circumstances warrant. This is perhaps the most 
elementary application of the seismic DCP concept. 

The following goals are assigned to the network: locate any earthquake within the Andean region 
of South America whose body wave magnitude is > 5.0 within two hours. The hypocenter must be 
located to ±0.1 degree in latitude and longitude and characterized as shallow, intermediate, or 
deep focus. 

The stations to be implemented were selected from the list of World Wide Standard Seismograph 
Network (WWSN) and array stations on the South American mainland. For the eight stations 
shown in Figure 15, a three minute P-S time circle will permit at least three stations to transmit P 
and S arrivals for events in the populous part of the western active area. Coverage is not so com* 
plete for the less populated areas of western South America (i.e., Tierra del Fuego) and the relatively 
inactive eastern area of Brazil. 

According to a study by Beirocal (1976), stations in continental South America observed 113 
events during 1973 with mb > 5.0 in the region bounded by latitudes 14®N and 56°S and longi- 
tudes 30°W and 90°W. Although fluctuations in this number occur on a yearly basis, the distrib- 
ution of stations and the quality of data used by Berrocal make it extremely unlikely that any 









events with mb > S.O were missed. Accordingly, the expected data load is about 90S station 
events per year. This will not tax the capability of any individual seismic DCP. 

To avoid transmitting unneeded events, each DCP would calculate a magnitude based on the usual 
log (A/T) (amplitude/period) calibrations to decide whether the event should be transmitted. 

Since the event buffer would haVe to be held during the calculation, the buffer swapping technique 
would be needed to avoid mssing an event. 

The GOES platform radio sets transmissions may be initiated in three ways; (1) The DCP may be 
polled by using interrogate channels; (2) the DCP may be activated by an internal timer on a regular 
basis, one or more times a day; and (3) the DCP may begin transmitting when a sensor threshold has 
been exceeded. These different modes are called; interrogate, self-timed, or emergency. 

After considering the objectives, the “emergency mode” was selected for the system baseline design. 
Advantages of the emergency mode over the other operation modes are; 

1. Requires a smaller DCP storage capacity. The emergency mode requires 108 kilobits, 
whereas the other modes require 324 kilobits for transmission on a 6-hour schedule. 

2. Requires a shorter playback time and dissipates lower power because the memory is 
smaller. 

3. Central station has near real-time monitoring of events. 

4. Requires only one master clock which is located at the central station. This is possible 
because the DCP delay can be measured during deployment; the transmission time through 
the GOES satellite system can be accountable, and the time from the event’s first zero 
crossing to event trigger pulse can be determined by the microprocessor and transmitted 
from the DCP with the event record. 
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The niiuor disadvantage of the emergency mode is that eight dedicated GOES channels are required 
versus one channel for the other modes. Also, an inoperative DCP could go unnoticed for several 
days. A combination mode, in which the E)CP returns housekeeping data once every six hours and 
operates in the emergency mode as well, may be the most desirable. 

The receive site requirements are modest. The DCS downlink from GOES is in the 1.7-GHz region. 
Microwave receiver technology for this kind of application is mature and the antennas are not very 
costly. Since the GOES spacecraft are in synchronous orbit, the ground antenna need only be 
positioned one time. Baseband signal processing is simple in the emergency mode since the exact 
baseband frequency for each DCP is known. A squelched discriminator can be used for each base- 
band signal with the squelch signal used to alert the data processing equipment. 

The data processing requirements are also modest. The basic analysis consists of two phases. In the 
first phase, the individual bit streams are converted to analog traces, timing information inserted 
and the traces are displayed for an analyst’s evaluation. At the same time, a preliminary hypocenter 
can be computed from the first arrival times reported by the DCP’s and the expected arrival times 
for other than main P can be marked. S arrivals, where present, would be selected. In the second 
phase, the analyst’s modifications would be used to calculate the final hypocenter and the individual 
traces would be output in final form. Neither of these tasks requires a particularly sophisticated or 
exp<>nsive computer. A microcomputer with video display, disk pack and hard copy plotter would 
be sufficient. Purchased now, the required hardware should cost much less than S20,000. 
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The DCPs are intended to operate unattended in the field. The only anticipated need for regular 
interventions would be the battery changes. This can not be avoided since the high power require- 
ments (10 W rf output for 18 minutes per buffer transmission) require high capacity batteries. Such 
batteries nonnally can not be recharged by solar panels. 

Field setup would include: 

a. Seismometer implacement and calibration. The DCP will need the necessary constants for 
a log (A/T) magnitude calculation. 

b. Processor activation and checkout. The noise characteristics of the site will determine 
the digital and analog constants. Since this will vary from site to site a certain amount of 
“cut and try” will be necessary. 

c. Transmitter activation, antenna pointing, and delay measurement. Antenna pointing 
angles can be calculated beforehand. The individual delay measurements can be made by 
initiating a transmission at a carefully measured time. The DCP clock, the delay transmis- 
sion initiation time, and the location of the E)CP (to a few meters) can be set by observa- 
tions of the Global Positioning System satellites during installation and activation. 

SUMMARY 

Our development effort has shown that there is no technical risk in building a Held worthy seismic 
DCP. Because of the advent of low power digital and analog electronics (in our case CMOS), a Held 
processor would require only a modest fraction of the total power budget. The power requirements 
of the DCP are dominated by the transmitter. 
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The nudor improvement over previous seismic signal processors is the use of the Allen (1978) 
event-recognition scheme. In the presence of severe cultural noise, Allen’s algorithm proved to be 
nearly perfect in its rejection of cultural signals. Our implementation on the CDPI802 microproc- 
essor required only modest amounts of ROM and RAM. 

The processor could be added to an existing GOES DCP with little difficulty. The digital format 
and bit rate costs nothing in terms of information content but does require a relatively long period 
of transmitter activation. The physical size of a field processor would most likely be about one 
quarter the size of an automobile battery. 

Operated in a dedicated network, through GOES, eight seismic DCPs would allow the location of 
large events (mb > 5.0) in South America within an hour of the event onset. Such a network 
would not tax the capabilities of the EXT* design and would represent no technical risk in its 
implementation. 
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APPENDIX A 

OPERATING INSTRUCTIONS 


Optrrttonal Proc»dur» 

The following describes the operational procedure for the seismic detector and the procedure to 
communicate with the system microprocessor via a silent 700-data terminal. 


1. Press power-on switch, located on power chassis, to apply power to seismic detector. 
Power lamp indicates power is applied to detector. 

2. Switch silent 700-terminal power on. 

3 . Switch helical-recorder power on. 

Rest 

Preu reset switch on seismic detector front panel. This places the microprocessor into initialization 
state and clears buffer memory. Event detect, record, and playback indicators are in off state. 

Run Utility Proarem 

1 . Press utility program switch. 

2. Depreu silent 700-termina] keyboard carriage return key, print head returns to left margin 
and types an asterisk. The asterisk acknowledges that the microprocessor address pointer 
is at memory location 0000. For programming instructions refer to “User Manual for the 
CDP1802 Cosmac Microprocessor,” RCA-MPM-201-B. 
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Run Srtimlc Program 


1 . Place unit into utility program state. 

2. Enter master program start address into memory location 0000 by typing keys 

iMOOOO COIOOO 

3. To verify, start address is stored into microprocessor memory type. 

?MOOO 3. 

4. Terminal response is to type 

0000 COlO 00. 

5. Press micro reset switch to enmre microprocessor is in initialization state. 

6. Press program start switch to initiate Rex Allen's seismic program. 

Seismic program will operate between two states (test and verify ) without additional co trolling. 
When an event is verified, the event indicator is turned “on" for approximately one second and the 
record indicator is turned “on" and remains on until memory buffer number one is filled. 

Buffer filling timing is predetermined by the status of switch E28 which is located on memory 
board (see Sheet 7). Filling time can be set from 10.2 to 1 53.6 seconds. (See Table A-1 for avail- 
able record times.) After the buffer is filled, the record indicator goes out and the playback indica- 
tor is turned on. Playback time is six times record time, if a second event occurs while buffer one 
is in playback state, the event is stored in buffer two and both the record and playback indicators 
are energized ; buffer two will play back after buffer one haii completed its playback. 
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Table A-1 

Buffer Memory Record and Playback 
Times and Switch Position 


Switch Position 

1 ime. Seconds 

Record 

Playback 

1 

10.24 

61.44 

2 

30.72 

184.32 

3 

51.2 

307.2 

4 

71.68 

430.08 

5 

92.16 

552.96 

6 

112.64 

675.84 

7 

133.12 

798.72 

8 

153.6 

921.6 


Program Constants 

Constants C| through are entered into program when microprocessor is in initialization state 
(reset). The programmed constants remain unchanged until the constant switch position is altered 
and the unit is reset. Table A-2 relates the constant switch positions to program value. 


Event Statistics 

Event peak values and time of occurrence referenced to first zero crossing are stored in the RAM 


memory. Memory location 0035 ‘‘ores the number of peaks. Memory location starting at 0100 
stores the time and peak values. There are three paired hexidecimal words per detected peak with 
the following format: 


3 paired hex words 


t 


(») 0 ) 0 ) 




P 

(n) (n) ^n) 
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Table A-2 
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11000 0.8 
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where n is the value recorded in location 0035, and i references i peak and runs from 0 to n. The 
event statistics can be obtained without disturbing the contents in the buffer. The procedures to be 
used to enter into utility run is made by using the micro reset switch. Events statistics have to be 
obtained from memory before the program goes into another test sequence. 

Holical Recorder 

Refer to Sprengnether VR-60 helical recorder operation manual for operation, calibration, and 
maintenance procedures. 

This is a three-channel recorder. Qiannel one records the real-time signal that has been amplified 
and filters to SO hertz. Channel two records the event-detect signal which is binary. Channel three 
records the delayed playback signal, which is 1/6 the rate of channel one signal. Channel three 
signal has been digitized, stored, and played back at the GOES channel rate of 1 00 bits per second. 

This signal is converted back to an analog signal before it is sent to the helical recorder. 

I 

j 


4 

1 
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APPENDIX B 

COMMENTARY ON THE CDP1802 
IMPLEMENTATION OF THE ALLEN ALGORITHM 

Event Detection Program Flow 

1 . Initialize and reset flags 
Read constants (C, • C, ) 
i= 1 

2. Input digital data R,, i = i 1 

( 1 ) convert to 2*s complement 

(2) R, = C,*R, 

(3) Calculations: 

R, = C,*(R,-R,.,) 

E, = R,» +AR,* 

3. Completed 2 second average? 

(i.e., i > 200) 

O.ltGo to 4 
0.2f Go to 2 

4. Compute reference level ( 7 j) 




5. Short term average abruptly increased? 

O.ltGo to 6 
0.2f Go to 2 

6. Save potential hit onset values 
0.1 T -i 

O 

0.2 = R, 

0.3D=R,-R,., 

0.4 M = 1 
0.5 S»0 

7. Save provisional peak 
P=|R,I 

8. Input digital data, i = i 1 

0.1 Convert to 2’s complement 
0.2Rj=C, *R, 

0.3 Calculations; 

R, = C, *(R,-R,.) 

E, = Rf +AR* 

9. Zero crossing? 

(i.e., R, = 0) 

O.ltGo to 11 
0.2fGoto 10 
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10. |R,|>P? 

O.ltToto? 

0.2fGo to 8 

11. 128 zero crossings recorded? 

(i.e.,M> 128) 

O.lt Go to 13 
0.2fGo to 12 

1 2. Record zero crossing 
01T„=i-T„ 

0.2A„=P 

0.3 M = M + 1 

13. Has 2 seconds passed since potential hit? 

(i.e.,i-T„>200) 

O.lt Go to 14 
0.2f Go to 8 

14. More than 40 zero crossings? 

(i.e., M > AO) 

O.lt Go to 15 
0.2fGoto2 

15. Declare significant event (set Q), and compute continuation criterion 
o = f(G,T^,M) 


16 . a,>o? 

O.lt S ■ 0, reset small count counter 
0.2f S » S + 1 

0.3 L B 4 i- M/4, value of S at which event is over 

1 7. Is the event over? 

(i.e.,S>L) 

O.lt Go to 18 
0.2fGoto8 

18. Declare event over (reset Q) 

0.1 i = 0 

0.2 Go to 2 
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APPENDIXC 

1802 ASSEMBLY CODE, MEMORY MAP, PROGRAM CONSTANT 
SPECIFICATIONS AND TIMING SUMMARY 
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0097 

• • 





1304 

; 

0098 

• * 
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1304 

8RI 

0088 

FHIFTi 8L0 

R8 


1305 

32219 

0100 

BZ 

EX1T8 

• • 

130? 

2fll 

0101 

BEC 

Rft 

■ • 

1308 

289 

0102 

DEC 

R8 

• • 

1309 

289 

0103 

DEC 88 



130^ 

F09 

0104 

LDX 


• • 

130B 

F69 

0105 

SHR 


• • 

130C 

A99 

0106 

PLG 

R8 


130D 

Ff^409 

010? 

RNI 

<H0 

.. GET PREV. SIGN 

130F 

32159 

0108 

BZ 

RHl 

.. SIGN ■ YESI 

1311 

889 

0108 

GLQ 

R8 

.. NDb IT"S 

1312 

F8809 

DUO 

□R1 

«S0 

.. RESTORE 

1314 

R89 

0111 

PLD 

R8 


1315 

889 

0112 

RHli GLO 

R8 

.. GET IT 

1316 

589 

0113 

STR 

R8 

.. S8VE HIGH ORDER BIT 

131? 

609 

0114 

IRX 



1318 

F09 

0115 

LDX 



1319 

?69 

0116 

SHRC 

* fl 

131f^ 

589 

Oil? 

STR 

R8 


131B 

609 

0118 

IRX 



131C 

F09 

0118 

LDX 



131D 

?69 

0120 

SHRC 



131E 

589 

0121 

STR R8 



131F 

30049 

0122 

BR 

RHIFT 

• • 

1321 

DFl 

0123 

EXIT8* SEP 

RF 

.. RETURN 

1322 

C012C09 

0124 

LBR 

HPLY 


1325 

9 

0125 

• • 



1325 

9 

0126 

• • 



1325 

• 

f 

012? 

• ■ 



1325 

9 

0128 

.. SUBROUTINE MEftT 

• • 

1325 

9 

0128 

.. CftLLINiS SEQ. •• 

• ■ 

1325 

9 

0130 


EX R4 . . 

R4« ADDCR) 

1325 

9 

0131 

• • 

• • 

R7* ft. oaEMR) 

1325 


0132 

SEP 

RE .. 


1325 


0153 

• m 



1325 

9 

0134 

m m 



1325 

9 

0135 

• • 



1325 

9 

0136 

LR»«06 



1325 

9 

013? 

Cl««08 



1 325 

9 

0138 

C2*«0C 



1325 


0138 

C3««0F 



1325 

9 

0140 

DELT8««25 



1325 

9 

0141 

CH8R-«28 



1325 

9 

0142 

TEMP2««3D 



1325 

9 

0143 

• • 



132*-; 

9 

0144 

• • 



1325 

9 

0145 

• m 



1325 

9 

0146 

□RG «1340 



1340 

F09 

014? 

HE8T: LDX 



1341 

5?9 

0148 

STR 

R? 

.. R7» R 

1342 

2?9 

0148 

DEC 

R? 


1343 

fe; 

0150 

SHL 


.. GET SIGN 

1344 

334FI9 

0151 

BDF 

ME688 

.. NEG YES 

1346 

F8009 

0152 

LDI 

0 

.. tiOf EXT ZEROES 

1348 

304C9 

0153 

BR 

ME8T8 

, , 

134A 

F8FF9 

U154 

NEG88: LDI 

XifF 

EXT ONES 

1 34C 

579 

0155 

ME8T8; SIR 

R7 

.. S8VE IT 

1 541* 

9 

0156 



.. TEMP= R f 16 BIT: 

1 3411 

1?9 

015? 

IHC 

R? 


134E 

F83D9 

0158 

LDI 

8. OaEMP2> 


1350 

889 

0158 

RLD 

R8 


1351 

F8085 

0160 

LDI 

Cl 


1 355 

889 

0161 

PLD 

R8 


1 354 

489 

0162 

LD8 

R8 



72 


13W 

BRI 

0163 

PHI 

RR 




09; 

0164 

LDh 

R9 


• 

1357 

rr; 

0165 

FLO 

RR 



1358 

£8; 

0166 

SEX 

R8 


• 

1359 

df; 

0167 

SEP 

RF 



135f=^ 

08; 

0168 

LBN 

R8 


• 

135B 

£41 

0169 

SEX 

R4 



135C 

731 

0170 

STXD 



. X« R4- ADD'LR) 

135D 

F800; 

0171 

LDI 

0 

• • 


135F 

R9; 

0172 

PLO 

R9 

• • 

SIGN DF LR 

1360 

rr; 

0173 

PLD 

RR 

• • 

SIGN OF R 

1361 

72; 

017^ 

LDXR 


• • 

A- LR 

1368 

FEI 

0175 

SHL 


ft • 

GET SIGN 

1363 

3368; 

0176 

BDF 

♦^5 

• • 


1365 

F8FFI 

0177 

LDI 

«FF 



1367 

R9; 

0178 

PLD 

R9 

• • 


1368 

Fo; 

0179 

LDX 


• • 

A« R 

1369 

fe; 

0180 

SHL 


• • 


136FI 

3B6F; 

0181 

BNF < 




136C 

f8Ff; 

0182 

LDI 

tiff 

• • 


136E 

rr; 

0183 

PLO 

PR 



136F 

24; 

0184 

DEC 

R4 

• • 


1370 

72; 

0185 

LDXR 


« • 

A« LR» H» R4» ADD'R' 

1371 

F5; 

0186 

SD 




1372 

57; 

0187 

SIR 

R7 

m • 


1373 

27; 

0188 

DEC 

R7 



1374 

89; 

0189 

GLD 

ri:9 

• • 

ADD EXTENSIONS 

1375 

57; 

0190 

STR 

R7 



1376 

E7; 

0191 

SEX 

R7 

• • 

P7=X= A. 1 aEMP;- 

1377 

8r; 

0192 

GLD 

RR 



1378 

74; 

0193 

RDC 


• • 


1379 

57; 

0194 

STR 

R7 


TEMP= R-LR» 16 BITS 

137R 

17; 

0195 

INC 

R7 



137B 

F825; 

0196 

LDI 

DELTR 



137D 

Rs; 

0197 

PLD 

R8 


• • 

137E 

F80c; 

0199 

LDI 

C2 
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1380 

891 

U199 

PLQ 

R9 



1381 

491 

0800 

LDA 

R9 



1388 

B8; 

0801 

PHI 

RH 

* * 

ftp. 1« LSC "C2" 

1383 

09> 

0808 

LDN 

R9 



1384 

881 

0803 

PLO 

RH 



1385 

E8I 

0804 

SEX 

R8 



1 388 

DFI 

0205 

SEP 

RF 

« • 

CALL MPLY<DELTA> 

1387 

041 

0806 

LOH 

R4 

• • 


1388 

FEI 

0807 

SHL 


• • 

SET SIGH 

1388 

338EI 

0808 

BOP 

NEGH 



138£ 

761 

0809 

SHRC 




138C 

30911 

0810 

£F POSP 



138E 

761 

0811 HEG8> 

SHRC 




138F 

FDOOl 

0818 

SDI 

0 



1391 

£81 

0813 PDS8t 

PHI 

RR 

• • 

Rfi.l- PBS*R> 

1398 

571 

0814 

SIR 

R7 

• • 

TEMP.O- PBS<R> 

1393 

F8891 

0815 

LDl 

CHRR 



1395 

881 

0816 

PLO 

R8 

• • 


1396 

F8001 

0817 

LDI 

0 



1398 

881 

0818 

PLO 

RR 

• « 


1399 

871 

0819 

DEC 

R7 



1398 

571 

0280 

STR 

R7 

■ • 


139B 

171 

0881 

INC 

R7 

a a 

R7« POO. OF MULTIPLICPNO 

139C 

DFI 

0888 

SEP 

RF 

• a 

CPLL MPLV < CHPR*RaR> 

139D 

F8841 

0883 

LDI 

R. OCDELTR- 

t> a 

• 

139F 

891 

0884 

PLO 

R9 

• • 


1380 

491 

0885 

LDA 

R9 



1381 

FEI 

0886 

SHL 


• • 


1388 

091 

0887 

LDN 

R9 



1383 

3B871 

0228 

BNF 

PDSB 

• • 


1385 

FDOOl 

0889 

SDI 

0 



1387 

£81 

0830 FDSB: 

PHI 

RR 



1388 

571 

0231 

STR 

R7 

• • 


1389 

871 

0232 

DEC 

R7 

• • 

T£MP« RBS<DELTR> 

1388 

F8001 

0833 

LDI 

0 



138C 

571 

0834 

STR 

R7 


• • 

138D 

171 

0835 

INC 

R7 



138E 

881 

0836 

PLO 

RR 


• • 

138F 

F83D« 

0237 

LDI 

DEMP2 



13B1 

881 

0238 

PLO 

R8 



13B8 

DFI 

0239 

SEP 

RF 

• • 

CPLL MPLV c TEMPcV* 

13B3 

5 

0240 .. 



* , 

= OELTPaac: 

13B3 

F8891 

0241 

LDI 

CHRR 



13B5 

891 

0242 

PLO 

R9 



1 3B6 

091 

0243 

LDN 

R9 



13B7 

F41 

0244 

HDD 




1 3B8 

591 

0245 

STR 

R9 



1 3B9 

891 

0246 

DEC 

R9 



13B8 

881 

0247 

DEC 

R8 
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13BB 

08> 

0848 

13BC 

741 

0848 

13BD 

581 

0850 

13B£ 

F889I 

0851 

13C0 

R8I 

0858 

13CI 

05 » 

0853 

i3ca 

F5I 

0854 

13C3 

57» 

0855 

13C4 

871 

0858 

13C5 

859 

0857 

i3ce 

889 

0858 

i3C7 

459 

0858 

13C8 

FBFF9 

088 C 

13Cft 

749 

088; 

13CB 

579 

0888 

13CC 

179 

0883 

13CD 

F80F9 

0884 

13CF 

R89 

0885 

13D0 

489 

0888 

13D1 

BR9 

0887 

13D8 

089 

0888 

1 3I»3 

RR9 

0888 

13D4 

F83D9 

0870 

13D8 

RSI 

0871 

13D7 

DF9 

0878 

13D8 

059 

0873 

13D9 

F49 

0874 

13DR 

559 

0875 

13DB 

859 

0876 

13DC 

889 

0877 

13IiD 

059 

0878 

13DE 

749 

0878 

!3DF 

559 

0880 

13E0 

159 

0881 

13E1 

DF9 

0888 

1 3E£ 

30409 

0883 

13E4 

9 

0884 

0000 




LDN 

R8 



RBC 




STR 

R8 

• • 

CHfiR- 

LDl 

CHRP 



FLO 

RS 



LDH 

R5 


. . «■ 

SB 


• 

. LO OFDEft 

STR 

R7 



DEC 

R7 


« • 

DEC 

R5 



DEC 

R8 



LDR 

R5 


• * 

KRl 

»FF 



RDC 



.. HIi3H ORDER BITS 

STR 

R7 



INC 

R7 


• m 

LDl 

C3 



RLD 

R8 


• • 

LDR 

R8 



PHI 

RR 


• • 

LDN 

R8 



PLO 

RR 


• • 

LDl 

TEMPS 



PLO 

RS 


• • 

SEP 

RF 

• • 

CRLl. MPLY 

LDN 

R5 



RDD 


• • 


STR 

R5 



DEC 

R5 

• • 


DEC 

R8 



LDN 

R5 

• • 


RDC 




STR 

R5 



INC 

R5 

• • 

R5» R.0‘RLFR> 

SEP 

RF 

• • 

RETURN 

BR NERT 




END 
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r 


MEMORY MAP 


4K 

RAM 


2.5 K 
PROM 


0.5 K 
ROM 


WORK SPACE 


ZERO CROSSINGS 
AND 

PEAK AMPLITUDES 
(MAX. 128) 


event detection 

TIMES 
(Max. 256) 



event detection 

PROGRAM 


M.l • 


-0100 


(HEX). 


• 1 ' •• 




-1000 




lUNUSED PROMS 




UT4 

UTILITY ROM 


REGISTER STORAGE 


-1400 

-1A00 

-8000 


-8200 

-8C00 

-8C=0 


OSto/N^r 
Of POOH 


f-'ci; 


IS 





High Pass Filter (Cl) 


I 

i 


Range; 0.5 to 1 .5 
Resolution: 0.1 
Mean: 0.995 

Steps: 1 1 


Index 

Value 

Real Value 

Index 

Value 

Real Value 

0 

0.5 

0.5000 

6 

1.1 

1.1093 

1 

0.6 

0.6250 

7 

1.2 

1.1875 

2 

0.7 

0.6875 

8 

1.3 

1.3125 

3 

0.8 

0.8125 

9 

1.4 

1.4062 

4 

0.9 

0.9060 

10 

1.5 

1.5000 

5 

1.0 

1.000 

11-15 

“NOT USED” 


Weighting Constant (C2) 

Range; 0.0 to 2.0 
Resolution: 0.15 
Mean: 0.65 

Steps: 14 


Index 

Value 

Real Value 

Index 

Value 

Real Value 

0 

0.00 

0.0000 

8 

1.05 

1.0625 

1 

0.15 

0.1406 

9 

1.20 

1.1875 

2 

0.30 

0.2812 

10 

1.35 

1.3750 

3 

0.45 

0.4375 

11 

1.50 

1.5000 

4 

0.60 

0.6250 

12 

1.65 

1.6250 

5 

0.75 

0.7500 

13 

1.80 

1.8125 

6 

0.90 

0.9060 

14 

1.95 

1.9375 

7* 

1.00 

1.0000 

15* 

2.00 

2.0000 


*Not required in original specification 
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1 


Range: 

0.2 to 0.8 


Resolution: 

0.1 



Mean: 

0.5 



Steps: 

7 



Index 

m 

Value 

Real Value 

0 


0.200 

0.2030 

1* 


0.225 

0.2187 

2* 


0.250 

0.2500 

3* 


0.275 

0.2730 

4 


0.300 

0.2960 

5* 


0.325 

0.3280 

6* 


0.350 

0.3437 

7* 


0.375 

0.3750 

8 


0.400 

0.4062 

9* 


0.425 

0.4218 

10* 


0.450 

0.4531 

11* 


0.475 

0.4687 

12 


0.500 

0.5000 


78 


Index Value Real Value 

13* 0.525 0.5234 

14* 0.550 0.5468 

15* 0.5750 0.5781 

16 0.6000 0.6093 

17* 0.6250 0.6250 

18* 0.6500 0.6562 

19* 0.6750 0.6718 

20 0.7000 0.7031 

21* 0.7250 0.7187 

22* 0.7500 0.7500 

23* 0.7750 0.7656 

24 0.8125 0.8125 

25-31 


“NOT USED” 


(C4) 


Range: 0.005 to 0.05 

Resolution: 0.005 
Mean: 0.025 

Steps: 10 


Index 

Value 

Real Value 

Index 

Value 

Real Value 

0 

0.0050 

0.0048 

10 

0.0300 

0.0312 

1* 

0.0075 

0.0078 

11* 

0.0325 

0.0332 

2 

0.0100 

0.0097 

12 

0.0350 

0.0351 

3* 

0.0125 

0.0126 

13* 

0.0375 

0.0371 

4 

0.0150 

0.0156 

14 

0.0400 

0.0390 

5* 

0.0175 

0.0175 

15* 

0.0425 

0.0429 

6 

0.0200 

0.0195 

16 

0.0450 

0.0449 

7* 

0.0225 

0.0234 

17* 

0.0475 

0.0468 

8 

0.0250 

0.0253 

18 

0.0500 

0.0507 

9* 

0.0275 

0.0273 

19-31 

“NOT USED” 




Threshold Constant (C5) 


Range: 

4.0 to 6.0 



Resolution: 

1.0 



Mean: 

5.0 



Steps: 

3 




Index 

Value 

Real Value 

0 

4.0 

4.0 

1 

5.0 

5.0 

2 

6.0 

6.0 

3 

7.0 

7.0 
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Expected Time Consumption 
(Search Mode) 



IR 1=127 


1250inst. = 10 ms) 

Program Seq. 

At Average 

Worst 

At Best 

or Function 

Values 

Case 

Values 

Germ Main 

130 

130 

130 

Sub Read 

20 

20 

20 

Sub Meat 

125 

125 

125 

Sub Mply (Cl) 

30 

320 


Sub Mply (C2) 

30 

245 


Sub Mply (R*) 

200 

200 

200 

Sub Mply (AR’ ) 

230 

230 

230 

Sub Mply (C3) 

50 

340 


Sub Mply (C4) 

300 

320 


Sub Mply (C5) 

90 

90 

90 

Totals 

1205 

2050 
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Expected Time Consumption 
(Validation Mode) 



|R|= 127 


(1250 inst. = 10 ms) 

Program Seq. 

At Average 

Worst 

At Best 

or Function 

Values 

Case 

Values 

Germ Main 

100 

100 

100 

Sub Read 

20 

20 

20 

Sub Meat 

125 

125 

125 

Sub Mply(Cl) 

30 

320 


Sub Mply (C2) 

30 

245 


Sub Mply(R*) 

200 

200 

200 

Sub Mply (AR*) 

230 

230 

230 

Sub Mply (C3) 

50 

340 


Totals 

815 

1660 
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